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By high temperature series expansion, exact diagonalisation and temperature density-matrix 
renormalisation the magnetic susceptibility x{T) and the specific heat C{T) of dimerised and 
frustrated S = 1/2 chains are computed. All three methods yield reliable results, in particular 
for not too small temperatures or not too small gaps. The series expansion results are provided 
in the form of polynomials allo-wing very fast and convenient fits in data analysis using algebraic 
programmes. We discuss the difficulty to extract more than t-wo coupling constants from the 
temperature dependence of x{T). 
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I. INTRODUCTION 

Quantum spin systems are amongst the most inter- 
esting and challenging problems of many-body theory in 
solid state physics. Due to their intrinsic many-body 
quantum character it is not possible to compute even 
simple quantities like magnetic susceptibilities or spe- 
cific heats in a straightfor-ward fashion. But there are 
by no-w a number of po-werful approaches like exact diag- 
onalisaton, quantum Monte Carlo, temperature density- 
matrix renormalisation or high temperature series expan- 
sion -which yield the desired quantities. 

The first aim of our present paper is to provide high 
temperature series data -which in combination -with some 
T = information can serve as an input for quick data 
analysis. High order series expansions constitute an ef- 
ficient, frequently used technique j^. Exact diagonali- 
sation and temperature density-matrix renormalisation 
-will serve as benchmarks to assess the reliability of the 
method proposed. Similar analyses are carried out in 
Refs. ^, ^ for unfrustrated dimerised spin chains and or- 
dinary spin ladders. 

The second aim of our work is to demonstrate that 
it is essentially impossible to deduce from one quantity 
like the magnetic susceptibility x(T) at not too small 
temperatures alone more than two of the three magnetic 
couplings of dimerised and frustrated spin chains. This 
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should caution anybody who is analysing such data in 
great detail. To illustrate the type of problem one can 
run in we refer the reader to the analysis of (VO)2P207 
which was analysed in the very beginning as dimerised 
chain Q . Then it was thought to be a two- leg spin-ladder 
with J| I w 1^ . But lately unambiguous evidence from 
inelastic neutron scattering was found jo that it is a set 
of weakly coupled dimerised chains ||^, The magnetic 
susceptibility x{T) is compatible with both scenarios 

In this article we focus on gapped one-dimensional 
spin systems. These form a rather large class com- 
prising dimerised spin chains, strongly frustrated spin 
chains but also spin ladders, cf. Fig. A representa- 
tive for a moderately dimerised spin chain is (VO)2P207 
^ H; a strongly dimerised spin chain realised in 



rised spin c 

Cu2(l,fC5Hi2N2)2Cl4 |lOl 0, |l| an example for 
a significantly frustrated spin chain is the spin-Peierls 
substance CuGeOa (see e.g. Ref. 14 and the discussion 
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therein) which is undimerised in its high temperature 
phase (T > Tc ~ 14K) but weakly dimerised in its low 
temperature phase. An important spin ladder compound 
is SrCu2 03 11. 

All these substances are interesting because they con- 
stitute disordered antiferromagnets with low coordina- 
tion number. This means that their ground state is 
not given by a Neel-type state (i.e. with finite sublat- 
tice magnetisation) but by a Resonating- Valence Bond 
(RVB) state made of superposed singlet-product states 
|l6| . If the systems are indeed gapped the average range 
of the singlet pairs present in the ground state is finite. 
In other words, the correlation length is finite. Gener- 
ally, an RVB state is favoured over a Neel state by low 
coordination numbers, by low values of the spin and by 
frustration (which simply weakens the classical ordered 
Neel-state). 

If it is possible to dope the insulating magnetic sys- 
tems unusual electronic properties emerge due to the 
strong interplay between charge and spin degrees of free- 
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dom. Some spin ladders like Sro.4Ca13.6Cu24O41.84 be- 
come even superconducting under pressure Of 
course, the appearance of a true phase transition at fi- 
nite temperatures requires a higher dimensionality than 
one 19, |20j. But the driving mechanisms can be 
present already in the low-dimcnsional systems. So, a 
deeper understanding of unusual electronic behaviour in 
doped antiferromagnets requires a thorough understand- 
ing of the magnetic subsystem. It is in this context that 
we perform the present investigation which is designed 
to determine the relevant magnetic couplings easily and 
reliably. 

Starting point of our theoretical study is the Hamilton 

N 

^ - E (-^ [(1 + (-1)''^) ^'^^+1 + "^^^^+2] ) (1) 

i=l 

with dimerised nearest and and uniform next-nearest 
neighbour interaction. The dimerisation is parameterised 
by 6. The ratio of nearest and next-nearest neighbor in- 
teraction is given by a. The dimerisation can arise from 
chemically different bonds as is the case in (VO)2P207. 
Alternatively, it may be induced by a static lattice dis- 
tortion via spin-phonon coupling as in CuGeOa. The 
Hamiltionian can also be viewed as a spin ladder with 
an extra diagonal coupling (1 — S)J (see Fig. p. In the 
limit (5 = 1 it is equivalent to a regular ladder model. 
In the limit S = 1, a = a system of isolated dimers 
is obtained. The ground state properties of the model 
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Figure 1: Dimerised and frustrated S = 1/2 spin chain. 
For S — 1 a. two-leg ladder is obtained. 



are investigated recently in numerous papers, see e.g. 

In the next section |^^e discuss briefly aspects of the 
methods we employ. The subsequent section II] is de- 
voted to the results for the magnetic susceptibil: 



lityx(7^)- 

In particular, we will discuss to which extent several cou- 
plings can be deduced reliably from x(T) data alone. 
Section |^ contains the results for the specific heat. We 
discuss how much can be learnt from C(T) data. The 
Summary ^ concludes this article. 



II. METHODS 

The present work is concerned with the finite temper- 
ature properties of the magnetic susceptibility x(^) ^^'^ 



the magnetic specific heat C{T). The analytic method 
we use is high temperature series expansion (HTSE). Its 
results are obtained as polynomials in the coupling pa- 
rameters with fractions of integers as coefficients so that 
no accuracy is lost. Details of the calculation are found in 
Ref. |2^. Furthermore, the data is provided in electronic 
form so that they can be put to use quickly. In order 
to maximise the range of applicability some extrapola- 
tion schemes are necessary which are described below in 
detail. 

In addition, we use numerical methods to cross-check 
the validity of the analytic results and to supplement 
them where necessary. A statistical method is quantum 
Monte Carlo (QMC) ||3^. In its present form, however, 
it suffers from the sign problem for frustrated systems 
so that we use it only in the unfrustrated cases. Mostly 
we employ the method of exact diagonalisation (ED) for 
smaller finite clusters It is based on the de- 

termination of all eigen- values so that the partition sum 
can be computed directly. Thereby any thermodynamic 
quantity is accessible. For not too low temperatures 
the results describe reliably the thermodynamic limit. 
The third technique that we use is temperature density- 
matrix renormalisation (T-DMRG)|^, By this nu- 
merical method very low temperatures can be reached 
reliably. 

We deal with the static magnetic susceptibility 



x(/3;5,a) = 
and the specific heat 
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In the Appendices the series coefficients are given. For 
the dimerised and frustrated chain the coefficients for 
both quantities x{T) and C{T) are provided up to order 
10 in the inverse temperature (3. For the unfrustrated 
dimerised chain they are given up to order 18. The bare 
truncated series, however, are not sufficient to describe 
the quantities under study at low values of T. To en- 
hance the region of validity of the high temperature series 
expansion Dlog-Pade representations [|4j are used. We 
show in the following that the complete susceptibility of 
gapped spin chains can be obtained if some low temper- 
ature information is incorporated in a simple way. This 
is explained subsequently in detail. 

A Dlog-Pade approximant of the expansion of the mag- 
netic susceptibility x{P) is given by 



(4) 



where Pi{[3') is the rational Dlog-Pade approximant with 
a polynomial of degree k in the numerator and a polyno- 
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mial of degree I in the denominator of 



where we use Pj' now for the rational function in u. 



We 



/ = dp, ln(4rx(/3')) = 



4Tx 



(5) 



Possible orders [k, I] of Piifi) have to fulfill fc + Z = n — 1 
where n is the order of the truncated series available. 
This can be seen by comparing the number of series co- 
efficients obtained and the number of independent (one 
less than the total number) coefficients in the two poly- 
nomials of Pj' {(3) . Due to the derivative on the right side 
of (||) one coefficient of the series expansion is lost. 

Additionally, we incorporate T = information to sup- 
plement the high temperature expansion. In the case 
of ungapped spin chains ||2^ we used knowledge of the 
full dispersion for this purpose. Though successful this 
approach is cumbersome in general. Therefore, we now 
simplify the procedure to incorporate T = information 
considerably. It is described for the case of gapped spin 
chains. 

At zero temperature the susceptibility of a gapped sys- 
tem vanishes and at finite but small temperature the de- 
viation is exponentially small due to the spin gap A 



X{T) 



for T < A . 
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Furthermore, the leading power in T can be determined 
on the basis of the dimensionality of the problem and of 
the behaviour of the dispersion close to its minima. For 
one-dimensional systems with quadratic minima , which 
is generic for gapped systems, one obtains |3q| 



for T < A 



(7) 



This equation provides information for two additional co- 
efficients of the Dlog-Pade approximant. Let us extend 
the series expansion of x by two terms a ■ /3"+^ -I- h ■ 
yielding Xc- Then the degree of the approximant Pi {(3) 
can be incremented by two fulfilling k -\- 1 = n -\- 1. The 
two additional conditions deriving from (|^) follow from 



-A - 1/(2/3) -fO(/3-2) 
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in the limit (3 oo. This condition requires that Pi {(3) 
is finite for (3 oo imposing the severe constraint I — k 
on the possible degrees of the Dlog-Pade approximant. 
To circumvent this constraint we substitute 



u = 13/(1 +(3)-^ 13^ u/(l - u) 



(9) 



thereby mapping the complete /3-interval [0, oo] to the u- 
interval [0, 1]. This mapping is justified by the continuity 
of Pl'{l3) in the limit (3 ^ oo. 

The asymptotic behaviour (|^) transforms under the 
mapping (0) to 



will henceforth consider only approximants in u so that 
no confusion should arise. 

The Dlog-Pade approximant P^ [u) is chosen such that 
it approximates du'^Txc/ ^Txc where (3 is expressed ac- 
cording to (^) as function of m. In this way, reliable in- 
terpolations between the low temperature behaviour and 
the high temperature series expansion can be obtained 
for arbitrary orders k and I complying with k + l = n -f 1 . 

To assess the range of validity of the Dlog-Pade ap- 
proximant various orders of Pj' are investigated. Exam- 
ination of the positions (3i of the poles of Pi and espe- 
cially of the largest modulus of them leads to a reliable 
estimate of the radius of convergence /9r of the series in 
j3 (Cauchy's theorem). We find roughly (3^ « 2/J as for 
the ungapped spin chains p9t . This means that the trun- 
cated series will always diverge around T w 0.5 J even for 
arbitrary high order in (3. Thus the polynomial represen- 
tation is not sufficient to describe the maximum of x{T) 
quantitatively. 



III. MAGNETIC SUSCEPTIBILITY x{T) 

A Applicability of the High Temperature Series 
Expansion 

In Fig. H various representations of the susceptibility 
are shown and compared. For fixed order I of Pj^ the 
Dlog-Pade representation of x moves for I = 2 upwards, 
for ^ = 4 downwards and for / = 5 in both directions (for 
I — i no evaluation is possible due to defective approx- 
imants in all orders). All representations converge for 
increasing order. The [5,4] and [4,5] approximants can- 
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Figure 2: Susceptibility for 5 = 0.1 and a — 0.24 for the 
dimerised, frustrated spin- 1/2 chain. Various orders of 
the Dlog-Pade approximant are shown. 
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not be discerned. But this is not a general feature. Other 
refiected approximants of the type [k, I] and [/, k] yield 
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Figure 3: Susceptibility of the dimerised chain for various 
values of S. The Dlog-Pade approximants are of order 
[15,4] 



O ED (N=16) 
--- T-DMRG 
HTSE 




T/J 



Figure 4: Comparison of high temperature series ex- 
pansion, exact diagonalisation and temperature density- 
matrix rcnormalisation data for a = 0.24 and S — 0.1. 



substantially differing results. In order to present our re- 
sults in a systematic and unbiased way we choose the P™ 
representation for all curves shown below (if not denoted 
otherwise) . Deduced from Fig. |^ we expect quantitatively 
reliable results down to T/J « 0.25 for the dimerised, 
frustrated chain. This conclusion is based on considering 
the highest orders and on looking for the range of T/J 
where they are consistent. 

In the purely dimerised case, see Fig. almost the 
whole temperature regime is excellently described. This 
is due to the high orders reached {0{/3^^)). In Fig. |[ the 
HTSE results are depicted in comparison to results from 
numerical methods (ED, QMC) and the exact result of 
the uniform chain pq |. In particular, the agreement be- 
tween the HTSE result and the exact one for the uniform 
chain is impressive. We think that this is the optimum 
which can be obtained by high temperature expansion 
since it is certainly not possible to assess the logarithmic 
low temperature corrections coming from the high tem- 
perature end. Technically, we used for the uniform chain 
instead of (fol) the obvious relations 
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The second relation ( |llh| ) reflects the fact that x(0) is 
finite. 

In Fig. I the [7,4] HTSE representation chosen is com- 
pared to exact diagonalisation and temperature density- 
matrix rcnormalisation data [^3| . The results are in very 
good accordance with each other. Only in the regime 
T/J < 0.2 there is a slight difference between the HTSE 
representation and the numerical results. 

In Fig. H the susceptibilities for various sets of param- 
eters are shown. The behaviour of the maximum of the 



susceptibility depends upon the parameters under study. 
For fixed next-nearest neighbour interaction a the posi- 
tion of the maximum moves to higher values of T/J for 
increasing i5 while the maximum value decreases. These 
effects are induced by the increasing gap. 

Fixing S, the position of the maximum moves to the 
left for increasing a. This can be understood from the re- 
duction of the dispersion on increasing frustration. The 
mobility of the excitations is more and more restricted 
[^ , [2^ , p8|| . The maximum value of x(T) remains almost 
constant. This can be seen as the result of an compen- 
sation of two contrary effects. On the one hand, the 
susceptibility would rise due to the shift of the maxi- 
mum position to lower temperatures where the global 
1/T factor (cf. (^)) enhances its value. But on the other 
hand, the frustration provides an additional antiferro- 
magnetic coupling in the system which works against an 
alignment of the spins. For instance, the antiferromag- 
netic next-nearest neighbour coupling induces a strong 
repulsion between aligned adjacent triplets on the dimers 



B Information Content of x(T) 

In this section we address the question to which ex- 
tent the parameters of the model Eq. |l] can be extracted 
from measurements of the susceptibility. In other words, 
we adopt the experimentalist's point of view who wants 
to determine the coupling parameters from experimen- 
tal data. Obviously, the main feature in the suscepti- 
bility curve is the maximum. So it is natural to use in 
the first place the maximum value Xmax and its position 
T = T,nax- We consider the product XmaxTmax since it is 
experimentally easily accessible and does not depend on 
the exchange coupling J. 
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Figure 5: Susceptibility for various values of S and a. 
Representation by the [7, 4] approximant except for the 
solid line in the uppermost panel where the [7, 4] approx- 
imant is defective and the [8, 3] approximant is used. 



The data shown in Fig. ^ are obtained by exact diago- 
nalisation of the full Hamitonian of a 16 site system. For 
a < 0.75 the data are exact to many digits for a system 
of this size. For a > 0.75 and (5 = finite-size effects of 
several percent occur. 

Note that for (5 = the quantity XmaxTmax reaches its 
minimum at a ~ 0.5 and then starts to increase again. 
This is due to the fact, that on growing a the system 
approaches two independent chains of half the size of the 
original chain with a = 0. Fig. |^ can be used easily: 
given the experimental input for XmaxTmax one can read 
off the value for 5 for a chosen a. 

To complete the analysis we plot in Fig. |^ the variable 
(/ /TiYi^-^ as a function of S for various values of a. So, once 
the value of S (for given a) is determined from Fig. 
Fig. helps to determine the exchange coupling J by 
reading off J/T^ax and multiplying by Tmax- 

We want to point out here that it is almost impossible 
to determine the values of J, a and d from the tempera- 
ture dependence of the susceptibility alone (cf. also 37). 
This phenomenon is well known from the investigations 
of (VO)2P207. In the case of this substance the suscep- 
tibility of the isotropic ladder i.e. a ladder with J|| = J±_ 
and the susceptibility of a dimerised spin chain with 
(5 = 0.2 fit both the experimental data equally well. For 
illustration the possible results for (VO)2P207, where 
Xmax — 2.07 10"'^ emu/mol V and T^ax — 74if corre- 



spond to the horizontal dashed line at 0.15 emu K/mol 
V in Fig. |6[ The difficulty to distinguish different sets 
of (J, (5, a) yielding the same value of Xmax^max is visu- 
alised strikingly in Fig. ^. The re-scaled susceptibilities 
belonging to various values of XmaxTmax are depicted. In 
the temperature region around the maxima and for larger 
temperatures the differences within each set are minute. 




Figure 6: XmaxT'inax versus 6 ior a = 

0.0,0.12,0.24,0.36,0.50,0.75,1.0,2.0 and 4.0 (in de- 
scending order at the right side of the graph, i.e. 
(5 = 1.0). The symbols are obtained by analysing the 
HTSE showing excellent agreement with the numerical 
data. For illustration, the dashed line refers to the 
experimental value of (VO)2P207. 



Thereby we conclude that it is impossible to determine 
all three coupling parameters from xiT) at moderate and 
at large values of temperature alone. In Fig. ^ the values 
of a and S are shown which belong to the various sets 
displayed in Fig. |[ 

Next we use the high temperature expansion results 
(see Appendix A2) to understand such a scaling be- 
haviour in the lowest orders of /3. It is obvious that the 
zeroth order 4rx ~ 1 does not allow the determination 
of any parameters. The first order of the susceptibility of 
a model with given frustration a and exchange coupling 
J is identical to the first order of another model with 
frustration ai and coupling Ji if 



Ji{l + ai) = J{l + a) 



(12) 



holds. In other words, if one had only first order results 
for the susceptibility it would be impossible to determine 
J and a independently. 

The second order of the HTSE depends on S. But 
again one can choose a particular value such that the 
sets (J, (2, (5) and (Ji, ai, Si) lead to identical zeroth, first 
and second order terms in /?. This choice is 



5l = 2ai ~ {2a - 



1 + ai 

1 + a 



(13) 



If the susceptibility is determined mainly by the first 
three orders the relations ( |l2| , p^ provide the recipe to 
re-scale the susceptibility such that different parameter 
sets yield similar temperature dependences. Indeed, if 
one focuses on the high temperature range this is true. 
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Figure 7: versus 5 for a = 

0.5,0.36,0.75,0.24,1.0,0.12,0.0,2.0 and 4.0 (in de- 
scending order at the left side i.e. (5 = 0.0). 



The precise position and value of the maxima, however, 
cannot be deduced from the first three orders alone. This 
implies that Eq. provides only rough estimates for the 
curves displayed in Fig. ^. 

To be cautious, we like to stress that we are not claim- 
ing that it is impossible to determine all three coupling 
J, a, 5 j/sufficient low temperature data is available to ob- 
tain the value of the gap. The gap A depends in a differ- 
ent way on the couplings than x(r) does |Q. Very often, 
however, the dependence at low values of the tempera- 
ture is either not accessible or it no longer corresponds 
to a pure ID system. In particular, interchain couplings 
lead to significantly altered gaps (for examples see 
Refs. |l^ although the behaviour at higher tempera- 
tures T > Jx is still well described by a ID model. 

IV. SPECIFIC HEAT C(T) 

It is a straightforward idea to extract further informa- 
tion about the magnetic properties of certain materials 
by considering also the specific heat C. Only in very rare 
cases, however, the magnetic part of C can be extracted 
in a reliable way from the measured data because the 
phononic contributions dominate C whenever the energy 
scale of the lattice vibrations is of the order of the mag- 
netic coupling J. 

At low temperatures where the phononic contributions 
vanish following the usual law a reliable extraction of 
C is possible. For high temperatures such a procedure 
will fail in general as can be seen, for instance, in a sim- 
ple model of Einstein (dispersionless) phonons coupled to 
Heisenberg chains |39|. So the full temperature depen- 
dence of the magnetic part of C can be measured only 
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Figure 8: Re-scaled susceptibilities for Xmax^inax^ 0.17, 
0.16, 0.15, 0.14, 0.13 emu K/mol for the a and 5 values 
shown in Fig. ^ (g-factor set to 2). 



in substances with a small exchange coupling J as it oc- 
curs, for instance, in organic magnetic materials, see e.g. 
Ref. 1^. 

Furthermore, it is in order to mention that there are 
indirect techniques to obtain the magnetic part of C(T) 
where the energy fluctuations are linked to dissipation. 
The latter is measured by the intensity of elastic scatter- 
ing in spectroscopic investigations, see e.g. The 
indirect approaches, however, may provide information 
on Tinax of C (T) but not on Cmax itself since overall fac- 
tors are not known. (In this section Tmax refers always to 
C{T). The position of the maximum of x{T) is denoted 




a 



Figure 9: Lines of constant height Xmax7max= 
0.17,0.16,0.15,0.14,0.13 emu K/mole in the a-5 plane. 
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In the Appendices A and B the coefficients for the spe- 
cific heat are provided. In order to compute C{T) addi- 
tional information on the low temperature behaviour is 
included as was done for x{T). At zero temperature the 
magnetic specific heat of a gapped system vanishes and 
at finite but small temperatures the deviation is exponen- 
tially small. For one-dimensional systems one has p55l 



C{T) oc 



(14) 



This asymptotic behaviour leads to a Dlog-Pade repre- 
sentation of C equivalent to the one of the susceptibility 
(see Eqs. | to |lC" 
reads 



The function approximated by 



9/3 In 



C(/3) 



(15) 



where the additional factor ensures that the argument 
of the logarithm tends to unity for /3 ^ 0. As for the 
susceptibility two additional terms are taken into account 
which are determined such that the correct asymptotics 
in u = /?/(l + /3) is guaranteed 



Pf{u)\ 
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(16) 
(17) 



problem arises since the maximum in the specific heat 
occurs at fairly low values of T/J « 0.3. In the purely 
dimerised case where we have reached 18*^ order in (3, 
the HTSE describes the maximum perfectly. Only in the 
regime T/J < 0.25 there is a slight difference to the ED 
and to the T-DMRG results. The same is true for the 
uniform chain (see inset of Fig. p^ ) where some devia- 
tions from the exact solution occur well below Tmax- We 
attribute these deviations below the maximum position 
2max to logarithmic corrections which cannot be assessed 
by the high temperature series expansion. 
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Figure 10: Specific heat C for various values of 6 and 
a. The Dlog-Pade approximants used are of degree [5, 4] 
for the frustrated and dimerised chains and [13, 4] for the 
chains without frustration. 



In Fig. IG the specific heat for various sets of parame- 
ters a anclT is compared to numerical ED and T-DMRG 
results. For large enough dimerisation 5, the ED and 
HTSE results agrees perfectly. For S — 0.1 and a — 0.24 
the HTSE describes the position of the maximum of the 
specific heat very well, but the absolute value deviates 
slightly from the ED and the T-DMRG results. This 



Figure 11: Re-scaled specific heat C for various values 
of S and a all yielding the same value XmaxT^ax = 0-15 
emu K/mol. 



Fig. pr] displays C(T) for the points in Fig. || belonging 
to XmaxT'max — 0-^5 cmu K/mol. Therefore, the temper- 
ature dependence is given in units of the maximum tem- 
perature TjJax of the susceptibility. Clearly, the curves 
differ from each other. Hence the knowledge of C(T), in 
addition to the knowledge of x(r), renders a complete 
determination of all three couplings possible. This is our 
main point in the present section. In other words, the 
knowledge of x(^) allows to fix 6 and J for given a. But 
a cannot be determined easily since there are sets of pa- 
rameters leading to very similar x(2^) curves, see Fig. |^. 
The corresponding C(T) curves, however, differ signifi- 
cantly as illustrated in Fig. ^ and thus provide a proper 
distinction of different parameter sets. 

To complete our analysis for the specific heat we pro- 
vide for C(r) in the Figs. |l| and ^ the analoga of Figs. | 
and for x{T)- Fig. |l^ displays the dimensionless (if fee 
is set to unity) specific heat which is independent of the 
value of the exchange coupling J. For given dimerisation 
d the frustration parameter a can be read off. Once a 
is known the curves in Fig. ^ allow to determine the 
energy scale. 

Let us briefly describe the principal behaviour of C(T) 
as a function of S and a. Some of the features of 
the curves can be understood by simple arguments. In 
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Figure 12: Maximum value Cmax of the specific heat as 
function of the frustration a for various values of the 
dimerisation i5. For i5 = no data is shown beyond a — 
0.5 since the exact diagonalisation results are not reliable 
there. 
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Figure 13: Position Tmax of the maximum of C{T) as 
function of the frustration a for various values of the 
dimerisation 5. For (5 = no data is shown beyond a = 
0.5 since the exact diagonalisation results are not reliable 
there. 



Fig. lU all curves converge to the dotted line for increas- 
ing values of a. The value of the dotted line is Cmax of a 
uniform chain without dimerisation and frustration. This 
is implied by the simple fact that the system approaches 
the limit of two independent chains for a ^ oo. Then 
the couplings J(l + 5) and J(l — 5) between the two legs 
(cf. Fig. |l|) become less and less important 



lim Cmax (a, 5) = Cmaxl" 0, (5 = 0) 



(18) 



For fixed value of the dimerisation 5 the position Tmax 
of the maximum of C is shifted to lower values on in- 
creasing a. This can be understood by the suppression 
of the dispersion of the elementary excitations due to the 
frustration |22, 2^, 28 . Thereby the overall energy scale 
on which excitations exist is reduced. The same was ob- 
served for the susceptibility as shown in the preceding 
section. A minimum is reached for a certain value of 
a ~ 0.5 because Tmax has to rise again since the system 
approaches the limit of two independent chains. Quanti- 
tatively, the large a limit fulfills 



lim 



Tlnaxla, 6) _ T'max(Q: = 0, 5) 



(19) 



Ja J 
independent of 5. 

It should be noted that for large 5 (cf. Fig. |l^) the 
value of Cmax does not substantially change which makes 
it difficult to discriminate the curves experimentally. The 
data shown are obtained for a TV = 16 cluster. A finite 
size analysis shows that the finite cluster results coincide 
with the infinite chain result except for points with (5 = 
and a > 0.5. 



V. SUMMARY 

The aim of the present article was two-fold. In the first 
place, we provided results and tools to facilitate and to 
expedite the analysis of experimental data in terms of a 
one-dimensional S* = 1/2 model, namely dimerised and 
frustrated spin chains. This model can also be seen as 
zig-zag chain and comprises in particular the usual spin 
ladder. Secondly, we demonstrated to which extent it 
is possible to determine the model parameters quantita- 
tively from the temperature dependences of the magnetic 
susceptibility x and of the specific heat C. 

We showed in detail how analytic high temperature 
series in high orders can be used to obtain reliable ap- 
proximants, namely Dlog-Pade approximants. The key 
point is to use additional well-known information on the 
T — Q and on the low-temperature behaviour to stabilise 
the approximants in the low-temperature region. We use 
the size of the gap, the form (linear or quadratic) of the 
dispersion in the vicinity of its minimum and the di- 
mensionality of the system as additional input. Thereby 
we achieve very good results in a straightforward fash- 
ion. The validity of the results is comparable to the one 
achieved by another extrapolation procedure introduced 
previously [ p9[ . The approach used in the present work is 
simpler since it requires less additional input. In Ref. |2^ 
knowledge of the whole dispersion was used. 

With the help of a computer algebra programme the 
approximants can be computed very quickly and easily. 
Thereby efficient data analysis becomes possible. 

The extrapolated series expansion results were gauged 
carefully by comparing them to numerical data. The 
methods employed are exact diagonalisation, quantum 
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Monte Carlo and temperature density-matrix renormali- 
sation. 

To ease data analysis further we included in thepresent 
work results for many sets of parameters. Figs. 0, |l^ 
and ^ make it possible to read off the coupling param- 
eters J, a and S if as little as the maximum values of 
the magnetic susceptibility, of the specific heat and the 
corresponding positions T^ax ^^'^ ^max known. 

It turned out that the knowledge of x{T) at moderate 
and high temperatures alone is not sufficient to deter- 
mine the three model parameters. Any additional knowl- 
edge, for instance on C{T) or on the singlet-triplet gap 
A, solves the problem. But such additional informa- 
tion is difficult to obtain. The specific heat is mostly 
dominated by the phonon contribution making it diffi- 
cult to be extracted. The gap A is in principle well de- 
fined. Frequently, however, the real systems lose their 
one-dimensionality at low energies, for instance due to 
small interchain couplings. Then the gap is influenced de- 
cisively by these additional residual couplings although 
the behaviour at moderate and higher temperatures is 
perfectly described by a one-dimensional model. In the 
analysis of experimental data it is certainly helpful to 
consider these facts. 
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APPENDIX A: DIMERISED AND FRUSTRATED CHAIN 
1 Specific heat 



Table I: Series coefficients 
frustrated chain C = ^ 


0'n,k,l 


for the high temperature expansion of the magnetic specific heat of the dimerised, 
a'^5'(/3J)". Only nonzero coefficients are presented. 


(n,k,l) 




(n,k,l) 




(n,k,l) 




(n,k,l) 


"n,k,l 


(n,k,l) 




(n,k,l) 




(2,0,0) 


3 
16 


(6,0,0) 


21 
4096 


(7,2,4) 


-819 
20480 


(8,4,0) 


35 
2048 


(9,3,6) 


1131 
36840 


(10,2,4) 


-276889 
5506024 


(2,0,2) 


3 


(6,0,2) 


-177 


(7,3,0) 


-413 


(8,4,2) 


-269 


(9,4,0) 


-261 


(10,2,6) 


51407 


16 


4096 


4096 


2048 


286720 


917504 


(2,2,0) 


3 
16 


(6,0,4) 


-333 
4096 


(7,3,2) 


661 
10240 


(8,4,4) 


1287 
10240 


(9,4,2) 


-111 
896 


(10,2,8) 


-214825 
3670016 


(3,0,0) 


3 


(6,0,6) 


73 


(7,3,4) 


1379 


(8,5,0) 


-407 


(9,4,4) 


21597 


(10,3,0) 


59305 


32 


4096 


20480 


61440 


286720 


2752512 


(3,0,2) 


9 


(6,1,0) 


63 


(7,4,0) 


651 


(8,5,2) 


2951 


(9,5,0) 


5901 


(10,3,2) 


38011 


32 


512 


20480 


61440 


81920 


393216 


(3,1,0) 


-9 
32 


(6,1,2) 


15 
128 


(7,4,2) 


637 
20480 


(8,6,0) 


-2449 
40960 


(9,5,2) 


-26673 
286720 


(10,3,4) 


-297061 
2762612 


(3,1,2) 


9 


(6,1,4) 


-123 


(7,5,0) 


-245 


(8,6,2) 


6943 


(9,5,4) 


537 


(10,3,6) 


-35787 


32 


512 


8192 


122880 


81920 


917504 


(3,3,0) 


3 


(6,2,0) 


-363 


(7,5,2) 


3297 


(8,8,0) 


1417 


(9,6,0) 


-2411 


(10,4,0) 


-138811 


32 


4096 


40960 


327680 


143360 


2752512 


(4,0,0) 


-15 


(6,2,2) 


303 


(7,7,0) 


917 


(9,0,0) 


-4303 


(9,6,2) 


-27 


(10,4,2) 


116115 


256 


2048 


40960 


688128 


143360 


1835008 


(4,0,2) 


3 


(6,2,4) 


-579 


(8,0,0) 


1417 


(9,0,2) 


-1401 


(9,7,0) 


-2229 


(10,4,4) 


5375 


128 


4096 


327680 


286720 


573440 


688128 


(4,0,4) 


-15 


(6,3,0) 


17 


(8,0,2) 


2199 


(9,0,4) 


3177 


(9,7,2) 


-11103 


(10,4,6) 


-125561 


256 


512 


81920 


573440 


573440 


1835008 


(4,1,0) 


-3 


(6,3,2) 


-45 


(8,0,4) 


5803 


(9,0,6) 


3531 


(9,9,0) 


-4303 


(10,5,0) 


51701 


32 


512 


163840 


286720 


688128 


1376266 


(4,1,2) 


3 


(6,4,0) 


105 


(8,0,6) 


11221 


(9,0,8) 


-18369 


(10,0,0) 


-334433 


(10,5,2) 


-14076 


32 


1024 


246760 


1146880 


110100480 


688128 


(4,2,0) 


-3 


(6,4,2) 


-69 


(8,0,8) 


-4997 


(9,1,0) 


2613 


(10,0,2) 


-249061 




-83087 


32 


1024 


983040 


573440 


22020096 


1376256 


(4,4,0) 


-IB 


(6,6,0) 


21 


(8,1,0) 


-4793 


(9,1,2) 


-7821 


(10,0,4) 


-133325 


(10,6,0) 


27641 


256 


4096 


61440 


143360 


11010048 


2752512 


(5,0,0) 


-15 
256 


(7,0,0) 


917 
40960 


(8,1,2) 


-1999 
20480 


(9,1,4) 


-5913 
286720 


(10,0,6) 


-33811 
3670016 


(10,6,2) 


3069 
57344 


(5,0,2) 


-15 


(7,0,2) 


1393 


(8,1,4) 


59 


(9,1,6) 


14211 


(10,0,8) 


-431449 


(10,6,4) 


-9087 


128 


40960 


4096 


143360 


22020096 


131072 


(5,0,4) 


-35 


(7,0,4) 


623 


(8,1,6) 


1981 


(9,1,8) 


-16347 


(10,0,10) 


49649 


(10,7,0) 


-1817 


256 


40960 


12288 


573440 


36700160 


917504 


(5,1,0) 


25 


(7,0,6) 


399 


(8,2,0) 


2323 


(9,2,0) 


3855 


(10,1,0) 


92629 


(10,7,2) 


-63813 


128 


8192 


24576 


57344 


2752512 


2762612 


(5,1,4) 


-25 


(7,1,0) 


-2611 


(8,2,2) 


-1347 


(9,2,2) 


-6261 


(10,1,2) 


6463 


(10,8,0) 


38993 


128 


40960 


40960 


286720 


172032 


1572864 


(5,2,0) 




(7,1,2) 


1281 


(8,2,4) 


-1369 


(9,2,4) 


-7923 


(10,1,4) 


-981 


(10,8,2) 


-324557 


T28 


40960 


40960 


57344 


458752 


11010048 


(5,2,2) 


-15 


(7,1,4) 


-2009 


(8,2,6) 


4679 


(9,2,6) 


40137 


(10,1,6) 


621 


(10,10,0) 


-334433 


128 


40960 


40960 


286720 


114688 


110100480 


(5,3,0) 


15 
128 


(7,1,6) 


3339 
40960 


(8,3,0) 


-59 
960 


(9,3,0) 


1 

10240 


(10,1,8) 


-206066 
2752512 






(5,3,2) 


-25 
128 


(7,2,0) 


-119 
4096 


(8,3,2) 


-41 
1280 


(9,3,2) 


1359 
35840 


(10,2,0) 


-420476 
11010048 






(5,5,0) 


-15 


(7,2,2) 


1491 


(8,3,4) 


671 


(9,3,4) 


-4311 


(10,2,2) 


6235 






256 


10240 


3840 


71680 


393216 
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2 Susceptibility 



Table II: Series coefficients an,k,i for the high temperature expansion of the magnetic susceptibility of the dimerised, 
frustrated chain x = ^ ; an,k,ioi''S''{(iJ)". Only nonzero coefficients are presented. 



t 1"! 


0'n,k,l 






^^I1,K,1^ 


0'n,k,l 


Cn If W 




i-n If W 

^n,K,i; 


0'n,k,l 




('•n,k,l 




1 


(0,3,0) 


-1 


(7,1,6) 


-3167 


(8,3,2) 


479 


(9,3,4) 


-74989 


(10,2,4) 


2661047 


4 


128 


1474660 


1290240 


27626120 


2972712960 


( ^ Ci fW 
(i,U,Uj 


-1 


(5,3,2) 


-1 


(7,2,0) 


-805 


(8,3,4) 


-3947 


(9,3,6) 


138421 


(10,2,6) 


-33937 


8 


384 


73728 


5160960 


82676360 


27525120 


(1,1,0) 


— 1 


(5,4,0) 


1 


(7,2,2) 


— 131 


(8,4,0) 


— 59 


(9,4,0) 


317 


(10,2,8) 


1754671 


~8~ 


512 


36864 


20160 


229376 


5945425920 


(2,0,2) 


— 1 


(5,5,0) 


— 7 


(7,2,4) 


— 857 


(8,4,2) 


5119 


(9,4,2) 


4099 


(10,3,0) 


— 311903 


Te 


5120 


368640 


2580480 


2949120 


82575360 


(2,i,Uj 


1 


(0,U,U ) 


-133 


(7,3,0) 


3023 


/' Q /I A \ 

(8,4,4 ) 


29 


/O A A\ 

(9,4,4) 


12337 


(10,3,2) 


-1295087 


8 


122880 


737280 


16360 


20643840 


743178240 


(3,0,0) 


1 


(6,0,2) 


— 83 


(7,3,2) 


1009 


(8,5,0) 


—877 


(9,5,0) 


—969 


(10,3,4) 


— 1311053 


96 


40960 


122880 


1290240 


655360 


743178240 


(3,1,0) 


1 


(6,0,4) 


— 21 


(7,3,4) 


— 583 


(8,5,2) 


—61 


(9,5,2) 


— 1387 


(10,3,6) 


943607 


T2S 


40960 


245760 


64512 


430080 


743178240 


(3,1,2) 


3 


(d,U,dJ 


-1129 


(7,4,0) 


-381 


{ Q a c\\ 
(o,D,Uj 


5389 


/n c a\ 

(9,5,4) 


69103 


/'in A n\ 

(10,4,0) 


9659 


128 


368640 


81920 


20643840 


41287680 


3932160 


(3,2,0) 


-1 


(6,1,0) 


9 


('1 A 0\ 

(7,4,2) 


113 


(8,6,2) 


6095 


(9,6,0) 


93463 


(10,4,2) 


3402433 


32 


1280 


147456 


4128768 


61931520 


1486356480 


( O O fW 

(^0,0, U ) 


1 


(D,l,2) 


1 


(7,5,0) 


943 


(8,7,0) 


-1271 


(9,6,2) 


-13629 


(10,4,4) 


-128473 


96 


1536 


368640 


1720320 


20643840 


41287680 


( A f\ (W 

(4,U,U) 


5 


(0,-1,4) 


3 


(7,5,2) 


-199 


\C).hX) ) 


i2G9 


(9,7,0) 


-67097 


(10,4,6) 


102007 


1536 


2560 


368640 


4587520 


82676360 


496462160 


(4,0,2) 


7 


(6,2,0) 


221 


(7,6,0) 


67 


(9,0,0) 


3737 


(9,7,2) 


5233 


(10,5,0) 


— 1177787 


768 


61440 


368640 


74317824 


9176040 


826763600 


(4,0,4) 


7 


(6,2,2) 


49 


(7,7,0) 


1 


(9,0,2) 


979 


(9,8,0) 


—361 


(10,5,2) 


52919 


612 


30720 


16128 


3440640 


1720320 


106168320 


(4,1,0) 


-23 


(a o a\ 

(6,2,4) 


-117 


(8,0,0) 


1269 


(9,0,4) 


481 


(9,9,0) 


3737 


/"in c a\ 

(10,5,4) 


153863 


768 


20480 


4587520 


1376256 


74317824 


495452160 


(A 1 o^ 


-3 


(0,0, U J 


-163 


(8,0,2) 


89 


(y,u,D ) 


263 


(10,0,0) 


-339691 


/" 1 n r\ \ 
(iU,D,U j 


599639 


256 


92160 


229376 


1474560 


5945425920 


594542592 


(4,2,0) 


1 


(6,3,2) 


97 


( 'H (\ A\ 

(8,0,4) 


1507 


(y,u,o ) 


15607 


(10,0,2) 


-215221 


/" 1 n o \ 
(10,6,2) 


-370969 


512 


30720 


20643840 


13762560 


5945425920 


212336640 


(4,2,2J 


7 


(6,4,0) 


7 


(8,0,6) 


-10831 


(9,1,0) 


-34337 


('inn A\ 
(10,0,4) 


195049 


/" 1 n A \ 
(10,6,4) 


-816989 


512 


15360 


10321920 


23592960 


2972712960 


2972712960 


(4,3,0j 


-1 


(6,4,2) 


-13 


(8,0,8) 


9623 


(9,1,2) 


-3899 


1 n n ci. \ 

(10,0,6) 


145961 


/" 1 n '7 n\ 

(10,7,0) 


791221 


96 


2560 


13762560 


1474560 


990904320 


1486356480 


(4,4,0j 


5 


(0,0, U ) 


23 


(8,1,0) 


-23629 


(9,1,4) 


-2423 


/'inn Q\ 
(10,0,8) 


1374211 


/'in'? o\ 
(10,7,2) 


14941 


1536 


7680 


20643840 


1835008 


1981808640 


99090432 




-7 


(0,D,Uj 


-133 


(o,l,2) 


36983 


(9,1,6) 


-12323 


/ 1 n n 1 n^ 

(10,0,10) 


-4776949 


/" 1 n o n\ 
(lU,o,Uj 


-367481 


5120 


122880 


20643840 


688 1280 


29727129600 


1486356480 


(5,0,2) 


1 


(7,0,0) 


1 


(8,1,4) 


76009 


(9,1,8) 


158933 


(10,1,0) 


-22843 


(10,8,2) 


-31027 


1536 


16128 


20643840 


165150720 


1486356480 


99090432 


(5,0,4) 


23 


(7,0,2) 


-59 


(8,1,6) 


-1927 


(9,2,0) 


14126 


(10,1,2) 


-322247 


(10,9,0) 


22433 


3072 


92160 


983040 


4128768 


247726080 


148636648 


(5,1,0) 


-49 


(7,0,4) 


-11 


(8,2,0) 


-58651 


(9,2,2) 


79 


(10,1,4) 


-16173 


(10,10,0) 


-339691 


6144 


9216 


13762560 


458752 


7741440 


5945425920 


(5,1,2) 


-67 


(7,0,6) 


-307 


(8,2,2) 


-4203 


(9,2,4) 


699 


(10,1,6) 


-14683 






3072 


92160 


656360 


458752 


9176040 






(5,1,4) 


7 


(7,1,0) 


5863 


(8,2,4) 


1061 


(9,2,6) 


22877 


(10,1,8) 


97039 






6144 


1474560 


393216 


20643840 


70778880 






(5,2,0) 


37 


(7,1,2) 


4289 


(8,2,6) 


33017 


(9,3,0) 


-1249 


(10,2,0) 


15205963 






1536 


491520 


41287680 


35389440 


5945425920 






(5,2,2) 


1 


(7,1,4) 


571 


(8,3,0) 


28751 


(9,3,2) 


-29209 


(10,2,2) 


515117 






512 


98304 


5160960 


9176040 


123863040 
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APPENDIX B: DIMERISED CHAIN 
1 Specific heat 



Table III: Series coefBcients an,i for the high temperature expansion of the magnetic specific heat of the dimerised 
chain C = ^„ ; UndS^pj)^ . Only nonzero coefficients for orders n > 10 in /3 are presented. 



(n,l) 




(n,l) 




(n,l) 


On, I 


(n,l) 




(n,l) 


On, I 


(11,0) 


37543 


(13,0) 


-1925339 


(14,12) 


-1410199439 


(16,6) 


-577431395917 


(17,14) 


28645978566427 


31457280 


41523609600 


543581798400 


846417258086400 


28566582460416000 


(11,2) 


-2943281 


(13,2) 


60223033 


(14,14) 


3603293663 


(16,8) 


-199750102373 


(17,16) 


-25549233744557 


1981808640 


39636172800 


41855798476800 


307788093849600 


228532659683328000 


(11,4) 


-1231703 


(13,4) 


2866603461 


(15,0) 


-31504270817 


(16,10) 


-141176079361 


(18,0) 


-80067486241427 


198180864 


871995801600 


362760263466600 


217660162079360 


4875363406577664000 


(11,6) 


-2446939 


(13,6) 


63260687 


(15,2) 


-7632646211 


(16,12) 


-6744074943121 


(18,2) 


5145779287709 


330301440 


19818086400 


10364292956160 


5078503548518400 


108341409036069200 


(11,8) 


-4199591 


(13,8) 


168008737 


(15,4) 


-89314894561 


(16,14) 


718080430229 


(18,4) 


8816140320217 


396361728 


68133063440 


72550050693120 


846417268086400 


36934671261962000 


(11,10) 


9792739 


(13,10) 


2484478763 


(15,6) 


-210219285347 


(16,16) 


-1288308081349 


(18,6) 


129606647760939 


1981808640 


435997900800 


217660162079360 


60942042682220800 


406280283881472000 


(12,0) 


3987607 


(13,12) 


-140182523 


(15,8) 


-3750596387 


(17,0) 


184265505341 


(18,8) 


211128422793049 


3170893824 


96888422400 


5580773130240 


3627502534656000 


812560567762944000 


(12,2) 


47342317 


(14,0) 


-369233453 


(15,10) 


-28368692533 


(17.2) 


501069785641 


(18,10) 


19613319318773 


13212057600 


930128855040 


51821464780800 


1904438830694400 


90284507529216000 


(12,4) 


70948027 


(14,2) 


-967322681 


(15,12) 


-110329916941 


(17,4) 


1750455145427 


(18,12) 


88360795513999 


26424115200 


1196879956480 


43530030415872 


5193924083712000 


406280283881472000 


(12,6) 


16985329 


(14,4) 


29175427 


(15,14) 


5917497137 


(17,6) 


65018005163 


(18,14) 


58,54118848751 


19818086400 


1993133260800 


1451()()l()i:iS(i24 


45:i4:!781fi8:i2f)f)f) 


8291434364928000 


(12,8) 


-238381 


(14,6) 


5006543507 


(16,0) 


85175S:i:i47()i 


(17,8) 


747l)!)42!)4U4!) 


(18,16) 


-859S()(i47(i98741 


1056964608 


5979399782400 


8706006083174400 


114266329841664000 


325024227105177600 


(12,10) 


2809583 


(14,8) 


2068244921 


(16,2) 


48804050567 


(17,10) 


-170193679963 


(18,18) 


6827887220393 


377487360 


1993133260800 


692523211161600 


5713316492083200 


1329644565430272000 


(12,12) 


-40555 


(14,10) 


1193462617 


(16,4) 


-448010792927 


(17,12) 


-221941384979 






117440512 


664377753600 


1171962357350400 


1269625887129600 
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2 Susceptibility 



Table IV: Series coefficients an,i for the high temperature expansion of the magnetic susceptibihty of the dimerised 
chain % = ^ J2n I '^n,lS\l3J)"' ■ Only nonzero coefficients for orders n > 10 in /? are presented. 



(n,l) 



(n,l) 



an,i 



(n,l) 



0-71,1 



(n,l) 



(n,l) 



(11,0) 
(11,2) 
(11,4) 

(11,6) 
(11,8) 
(11,10) 



54499737600 

-2683 
29030400 

-674059 
7431782400 

-56827 
1238630400 

98731 
707788800 

-9539 



(12,0) 
(12,2) 
(12,4) 
(12,6) 
(12,8) 
(12,10) 
(12,12) 



2242274918400 

-35760853 
2615987404800 

-50519537 
1046394961920 

-18989863 
356725566200 

-212653873 
6231974809600 

-161397977 
623197480960 

193615997 
5231974809600 



(13,0) 
(13,2) 

(13,4) 
(13,6) 
(13,8) 
(13,10) 
(13,12) 



809710387200 

2790083 
1189()85184()() 

9861031 
653996851200 

16757 
19818086400 

-56229473 
3923981107200 

-9467111 
71345111040 

2172449 

21799895040 



(14,12) 
(14,14) 



(14,0) 
(14,2) 
(14,4) 
(14,6) 
(14,8) 
(14,10) 



3957275492352 

361148659 
35876398694400 

9492225643 
466393183027200 

5394474767 
296795661926400 

7897101007 
652950456238080 

-10743048697 
652950456238080 



(15,0) 
(15,2) 
(15,4) 
(15,6) 
(15,8) 
(15,10) 
(15,12) 
(15,14) 



378600476623 
3264752281190400 

-656281799 
76947023462400 



-65174099663 
28566582460416000 

-49949749 
11719623573504 

373889611 
357082280755200 

196672854197 
34279898952499200 

1134039913 

173130802790400 

86433646963 
11426632984166400 

125245014907 
1713994947624960 

-104930723893 
3808877661388800 



(16,0) 
(16,2) 
(16,4) 



-258645079463 
498616712036352000 

-15175092143 
3767021862912000 

-8706451935593 
1371 195958099968000 



(16,6) 
(16,8) 

(16,10) 
(16,12) 
(16,14) 
(16,16) 



-381013820701 
85699747381248000 

-6536122273267 
2742391916199936000 

176252096491 
342798989524992000 

32635888627387 
1371195958099968000 

-150602416817 
3808877661388800 

10785364513867 
5484783832399872000 



(17,14) 
(17,16) 



(17,0) 
(17,2) 
(17,4) 
(17,6) 
(17,8) 
(17,10) 
(17,12) 



1228965600979 

2590036809744384000 

146506668199 
685597979049984000 

-134460640739 
62327089004544000 

-804378832603 
228532659683328000 

-1061886885637 
342798989524992000 

-535584743809 
228532659683328000 

344293742347 
228532659683328000 



-22234161829843 
685597979049984000 

52097662147 
7031774144102400 



(18,0) 
(18,2) 
(18,4) 
(18,6) 
(18,8) 
(18,10) 
(18,12) 
(18,14) 
(18,16) 
(18,18) 



1116582102301823 
4475583607238295552000 

33268250832001 
27124749134777548800 

301983779893 
208944146996185600 

30860378761391 
74693060120638259200 

-182357607317269 

745930601206382592000 

-114706799537333 
248643533735460864000 

-<)()44275()894411 
74593060120638259200 

-169097404668739 
10666161446806465600 

6285092004168191 
497287067470921728000 

-468309667465837 
1032826986285760512000 



